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SUMMARY 

Delamination growth in compressively loaded composite laminates was 
studied analytically and experimentally. The configuration used in the study 
was a laminate with an across-the-width delamination. An approximate super- 
position stress analysis was developed to quantify the effects of various geo- 
metric, material, and load parameters on mode I and mode II strain energy 
release rates Gj and G.^, respectively. Calculated values of Gj and G were 
then compared with measured cyclic delamination growth rates to determine the 
relative importance of Gj and G^. High growth rates were observed only when 
G I was large. However, slow growth was observed even when G was negligibly 
small. This growth apparently was due to a large value of G . 

INTRODUCTION 

In composite structures subjected to compression loads, delaminations can 
cause localized buckling (fig. 1). High interlaminar stresses at the edges of 
the buckled region often lead to cyclic delamination growth (herein referred to 
as instability-related delamination growth) . 

The objective of this paper is to investigate the mechanism of instability- 
related delamination growth. Figure 1 shows the configuration used in the 
study a laminate with a "through-width" delamination. This configuration 
was selected because it is perhaps the simplest configuration that exhibits 
instability-related delamination growth. Goals of the investigation were: 

(1) to develop and use an approximate superposition stress analysis to explain 
how various geometric, material, and load parameters affect interlaminar 
stresses, (2) to determine the delamination growth behaviors predicted by 
several different criteria based on strain energy release rates, and (3) to 



compare analytical calculations with experimental observations to determine 
the applicability of each growth criterion. 

Because of the stress singularity at the end of the delamination (crack 
tip), calculated stresses there have little meaning. Strain energy release 
rates are finite parameters which characterize the intensity of the stresses 
near the crack tip. Consequently, in the following discussion strain-energy 

release rates will be used to characterize the severity of the interlaminar 
stresses. 
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NOMENCLATURE 

half-length of delamination before loading 
half-length of delamination after loading 
virtual crack closure distance used in strain energy 
release rate calculations 
specimen width 

arbitrary constants 

bending stiffness of the buckled region given by 
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where p = number of plies 


d x> d y unit load solutions for displacements near crack tip 

k 

E Young's modulus for ply k 

E ll> E 22> E 33 Young's moduli of unidirectional ply. The subscripts 1, 

2, and 3 refer to the longitudinal, transverse, and 
thickness directions respectively. 
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unit load solutions for forces at crack tip 
Mode I strain-energy release rate 
Mode II strain-energy release rate 

maximum possible value of Gj for current delamination 
length 

shear moduli of unidirectional ply 
moment 

number of applied load cycles 

axial loads in regions A, B, C, and D respectively 
remote applied compressive load 

axial stiffness of regions A, B, C and D given by 
P 

S = b£ E k (a - a ) 
k=l k k_1 

where p = number of plies 

thickness of buckled region 

rectangular Cartesian coordinates 

distance from top surface of laminate to ply "k"; 

top ply is ply 1 

lateral deflection at x = -a due to applied load 

/N 

value of 6 corresponding to Gj 
initial lateral deflection at x = -a 
Poisson's ratios for unidirectional ply 
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ANALYSIS 


The configuration shown in figure 1 was idealized as a two dimensional 
plane strain problem. Linear and nonlinear finite element analyses and an 
approximate superposition analysis were used to calculate strain energy 
release rates for the two dimensional idealization. The nonlinear analysis 
was used to provide reference solutions for evaluation of the approximate 
superposition analysis . The linear analysis was used to calculate several 
constants used in the approximate superposition analysis. The nonlinear 
analysis is described in reference 2, and the linear analysis is simply a 
linear version of this analysis. 

The approximate superposition analysis, the procedure for calculating 
strain energy release rates, the finite element models, and material 
properties are discussed in the following sections. 

Approximate Superposition Analysis 

Superposition techniques have been widely used in linear stress analysis 
to represent a complicated problem as a combination of several simpler problems. 
Application of the principle of superposition to nonlinear problems first 
requires a transformation that results in a linear system. 

The key to the transformation is replacement of the source of nonlinearity 
with equivalent loads (fig. 2a and b) . Because of symmetry only half of the 
configuration is considered. The buckled region (which responds nonlinearly 
due to significant rotations) is replaced by the loads P D and M, the axial 
load and moment respectively in the column where it is cut (fig. 2b). The 
new configuration is linear, with three nonlinearly related applied loads 

P T» P D > and M * B y superposition the number of loads can be reduced to two, as 
illustrated in figures 2c - 2e. 
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The load system in figure 2c is divided into the two load systems shown 
in figures 2d and 2e. Because is calculated using rule of mixtures, 

the load system in figure 2e causes a uniform axial strain state and no inter- 
laminar stresses. Consequently, in terms of interlaminar stresses, only the 
load system in figure 2d (ie. (P^) and M) need be considered. Accordingly, 
in the current study involving strain energy release rates, figure 2d is 
the linearized equivalent of the nonlinear problem in figure 2a. 

The appendix describes a strength of materials analysis for calculating 
< P C “ P D^ anc * T ^ e key equations from the appendix are 


( 1 ) 

( 2 ) 

(3) 

To use the loads (P c - Pp) and M in a two dimensional analysis requires that 
they be expressed as an equivalent distribution of tractions. To calculate 
this distribution, the axial strains were assumed to vary linearly through the 
thickness where the tractions are applied (ie. at the cut). Intuitively, this 
seems to be reasonable if region D (fig. 2) is not cut too close to the 

crack tip. The validity of the assumed linear variation will be checked later 
in this paper. 

Linear finite element analysis was used to calculate the response of the 

linearized configuration in figure 2d to unit values of (P - P ) and M. 

CD 

Because the configuration is linear, the solution for any arbitrary combination 
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of (P c - P D ) and M is simply a linear combination of the unit load 
responses. If region B (fig. 2) is much thicker than region C, the unit 
load solutions are very insensitive to delamination length. In the current 
study the ratio of thicknesses was 61 to 3. Hence, the unit load solutions 
for 2a = 25 mm were used for analysing all delamination lengths. Also initial 
waviness of the buckled region does not enter into the finite element analysis. 
Delamination length and initial waviness were both accounted for in the 
strength of materials analysis in calculating (P c - P D ) and M, equations 

(2) and (3) respectively. This procedure will be discussed further in the 
next section. 

Strain Energy Release Rate 

The virtual crack closure method (ref. 1) was used to calculate mode I 
and mode II strain-energy release rates, Gj. and respectively. The forces 
transmitted through the node at the crack tip and the relative displacements 
of the two nodes on the crack boundary closest to the crack tip were used in 
the calculation. Equations 4 show how this technique is used for the super- 
position stress analysis. 

G I = 2Zab < P C-V F y +MI >j [< P C ' V d J + Md y 

J J (4) 

G II - 2S3b [ < p c - Vx + m l ] [< P C - V d * + M X 

In these equations F x , Fy, d x , and dy are the unit load values of the nodal 
forces and the corresponding relative nodal displacements in the x and y 
directions. (The coordinate system is defined in fig. 2.) The superscripts 
1 and 2 on the unit load parameters identify parameters associated with 
“ Pp) and M, respectively. 
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If the distance is small between the crack tip and the nodes used to 

calculate relative displacements, then F^/d^ ** F^/d^ and F^/d^ = F^/d^ 

y y y' y r x /a x x /a x‘ 

Using these relationships in eqns. (4) results in 

y 1 J 

d i (5) 

G n = ihb [ (p c ~ V F i + m l\ 

x L J 

In the results and discussion section it will be shown that for high loads 
or long delamination lengths, Gj. is zero, i.e., the crack tip closes in the 
normal direction. To prevent the crack faces from overlapping (analytically) 
requires the addition of multipoint constraints on the crack face nodes. Con- 
ceptually, the crack face nodes are connected in the direction normal to the 
crack face by infinitesimal springs. These springs have infinite stiffness 
in compression and zero stiffness in tension. To determine whether to select zero 
or infinite stiffness requires solution of a nonlinear contact problem. To 
include the contact problem directly in the superposition analysis would 
severely complicate the otherwise simple equations. Therefore, use of a non- 
contact analysis to approximate G n was investigated. 

A laminate with 2a = 76.2 mm was analyzed using two different approaches. 
First contact forces were ignored (i.e., overlap of crack faces was allowed). 

Gj. and G ri were calculated using equation (5) . In the second approach, overlap 
of the crack faces was prevented, which is more realistic. G n was calculated 
using equation (5). (Note that eq. (5) yeilds Gj = 0 when overlap is prevented.) 
Applied loads (P T ) ranged from 14.8 kN, which corresponds approximately to 
initial crack tip closure, to 55.2 kN. 
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When crack face overlap was prevented, a larger value of Gjj- 
calculated than when overlap was allowed. The difference in the Gjj values 

increased with load. But in all cases the difference was approximately equal 
Gj calculated using the approach which allowed crack face overlap. For 
example, for P T » 55.2 kN the contact analysis yielded G TI = 413 J/m 2 . When 
crack face overlap was allowed, Gj. and G Ir were 35 and 384 J/m 2 , respectively. 
The sum of these values is within approximately 1.5 percent of the more 
realistic solution, i.e., G^j = 413 J/m 2 . Apparently the crack-face contact 
forces do not significantly alter the total strain energy release rate. Hence, 
when there is crack tip closure, the total strain energy release rate from the 
non-contact analysis can be used to approximate G XI (which is then the total 
strain energy release rate, since Gj is identically zero). 

Finite Element Model 

A typical finite element mesh for the nonlinear analysis is shown in 
figure 3. Because of symmetry only half of the laminate was modeled. The 
mesh contains 813 nodes and 740 four-node isoparametric elements. Reduced 
integration was used to improve the performance of the elements in modeling 
bending deformations. Because the rotations are small except in part of .the 
buckled region, the nonlinear strain— displacement relations were used only for 
the region y > 0, -a < x < -0.56 mm. As shown in figure 2(b), the 
linearized configuration is the same as the nonlinear configuration except that 
most of the buckled region is removed. Accordingly, the mesh used in the 
linear analysis was derived from that in figure 3 by removing elements in the 
deleted part of the buckled region. 
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Materials Properties 

The material studied was NARMCO T300/5208* graphite/epoxy. The 
unidirectional ply properties were assumed to be 

E 11 = 140 GPa 
E 22 = E 33 = 14 GPa 

v 12 - v 13 * v 23 ’ °- 21 
°12 ' g 13 * G 23 ' 5 ' 9 GPa 

Plane strain (i.e,, = 0) and E “0 were imposed to calculate the 2D 

properties. In regions where coarse finite elements spanned several plies, 
laminate theory was used to obtain average properties. 


*Use of trade names or manufacturers does not constitute an official 
endorsement, either expressed or implied, by the National Aeronautics and 
Space Administration. 
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EXPERIMENTAL PROCEDURE 


The specimens used for this study were fabricated and tested by Northrop 
Corporation. (Details will appear in NASA CR-166046, "Performance of a Quan- 
titative Study of Instability-Related Delamination Growth," by R. L. Ramkumar.) 

A cursory description of the experimental procedure is given herein. 

The specimen consisted of 64 plies of T300/5208. The fiber orientation and 
stacking sequence were r0 4 / (0/45/90/-45) ? J g . The laminate width, b, was 25.4 mm 
To simulate a delamination, kapton film was used to prevent bonding over a 
19 mm length between the third and fourth plies. The ply thickness was 
assumed to be 0.14 mm. Six specimens were tested in fatigue under 
compressive constant-amplitude loads. Minimum compressive load was ten per- 
cent of the maximum compressive load. The load frequency was 10 Hz. Delami- 
nation lengths were measured with a microscope. 

RESULTS AND DISCUSSION 

First, the accuracy of the approximate superposition analysis will be 
evaluated. Then the effect of various parameters on Gj. and Gjj will be con- 
sidered. Finally, the experimental observations will be compared with the 
analytical results. 

Evaluation of Approximate Superposition Analysis 

The approximate superposition analysis was evaluated by comparison with 
results from a geometrically nonlinear finite element analysis. Recall that a 
major assumption in the approximate analysis was that the strains vary linearly 
through the thickness where (P c - P D ) and M are applied. Fig. 4 shows the 
axial strain variation through the thickness at x « - 0.2 and - 0.7 mm obtained 
using the nonlinear finite element analysis. Along the line x = -0.7 mm the 
strains vary almost linearly for the three applied loads. However, closer to 
the crack tip along x = - 0.2 mm the variation is more nonlinear, especially 
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near y - 0. In the following, (P c - P D ) and M were applied at 
x = -0. 76 mm. 

The unit load solutions F F d. and d are 

* y x y 

F* « 9.36 x io“ 5 

X 

F* « Q. 531 nf 1 

fJ - Q.0261 

Fy = -0.252 nf 1 

d* = 1.40 x io -10 m _ N -1 

dy = 1.17 x 10" 8 N _1 

These values were used in equation 5 for any combination of (P - P ) and M 

C D 

to obtain Gj and G--. 

Differences between Gj and Gjj from the approximate superposition 

analysis and the geometrically nonlinear finite element analysis can be traced 

mainly to two sources: (1) nonlinear variation of the strains through the 

thickness of the buckled region and (2) inaccuracy in determining (P - P ) 

C D 

and M. 

By using (P c - P D ) and M from the geometrically nonlinear finite 
element analysis, the effect of nonlinear variation of the strains can be 
examined. Figure 5 shows that this effect is small. 

Figure 6 shows that if the strength of materials analysis is used to calcu 
late (P^, - Pp) and M, the difference is much larger. Hence, most of the 
difference between the two analyses is due to inaccuracy in determining 
( P C “ P D^ and M * But the Seneral trends for the Gj variation with 
delamination length and load are predicted very well. In figure 6 the curves 
for the two analyses seem to differ (approximately) by a constant scale factor. 
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A direct test of the approximate analysis for predicting trends is to use 
it to coalesce the curves in figure 6 into a single curve. Equations (2) and 
(3) show that (P^ - P^) and M can be expressed as functions of 6. Hence, 

from equations (5), and are functions of 6. Equations (1), (2), 

_2 

(3), and (5) show that for constant 6, P T varies as a and G^ and G^ 
-4 4 2 

vary as a . Hence, plotting a Gj. vs a P^, should coalesce the curves for 
various delamination lengths. Figure 7 shows that the data for five delami- 
nation lengths (including those in fig. 6) do coalesce into a narrow band 
around a single curve. Since the peak values of Gj. for various lengths 
differ by more than two orders of magnitude, the closeness of the fit suggests 
the approximate analysis is accurate for predicting trends. Therefore, all 
results that follow are obtained with the approximate superposition analysis. 
Figure 7 also shows that if nonlinear finite element results are available for 
one delamination length, the values for other lengths can be estimated 
immediately . 

An advantage of the superposition analysis is that it allows a problem to 
be dissected. In particular, one can determine the relative importance of the 
loads (P^ - Pq) and M on G^ and Gj.j. Figures 8 and 9 show G^. and G^j 
calculated by using M alone and by using (P^ - P^) and M in combination. 
Although intuition might suggest that only the peeling action caused by the 
moment M has a significant effect on G^, figure 8 shows that (P^ - P^) 
contributions cannot be ignored. Figure 9 shows that both (P c - P^) and M 
are also important when calculating Gjj. 

Parametric Study 

The effects of several parameters on Gj. and G^ were examined using 
the approximate superposition analysis. The parameters were initial waviness, 
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delamination length, applied load, and the ratio of axial to bending stiffness 
for the buckled region. 

Initial imperfections in the form of simple sinusoidal waviness were 
assumed (eq. (6)) 


v(x) 

INITIAL 



( 6 ) 


where v(x) = Distortion in the y direction. When a column is initially 
wavy, bifurcation buckling does not occur. As soon as load is applied, the 
column begins to deflect laterally, which causes interlaminar stresses. Hence, 
Gj and Gjj are nonzero as soon as load is applied. If 6 q = 0, Gj and G IJ; 
are zero until buckling occurs. However, figure 10 shows that the peak value 
of Gj is significantly reduced, even for very small imperfections. In 
contrast, figure 11 shows that G^ is hardly affected by initial waviness. 

Figures 6, 12, and 13 show the effect of delamination length on G^ and 
<»Ij. The shorter delaminations have the larger values of peak Gj. (i.e., Gj.) 
(Us* 6). However, for the longer delaminations G^ becomes nonzero at lower 
loads. Figure 12 shows that after only a little delamination growth, Gj 
reaches a peak and decreases rapidly with further growth. At 2a = 40-50 mm, 
the crack tip closes in the normal direction and Gj is identically zero. 
Further delamination growth causes compressive normal stresses to develop at 
the crack tip. In contrast, G^^ initially increases then decreases only 
slightly to a constant value with increased delamination length (fig. 13). 

Note that G IX is typically much larger than Gj. 

Figures 6 and 14 illustrate the effects of applied load on Gj and 
Gjj, respectively. The mode I strain energy release rate G^ first 

/s 

increases to a peak value (G^), then decreases with increasing load 
(fig. 6). In contrast, G^j monotonically increases with increasing 
load. As a result, G^ and Gj^ do not usually reach peak values 
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at the same time during a fatigue load cycle. Furthermore, the load at which 
Gj is maximum decreases with increasing delamination length. 

A 

Since Gj is the maximum possible value of Gj for a given delamination 
length, it is of interest how Gj varies with delamination length. The first 
step in determining this variation is to determine the lateral deflection <5 

Ss A 

corresponding to G^.. The lateral deflection 6 is obtained by solving 



(7) 


Equations (2), (3), (5), and (7) are combined to obtain the governing 
equation (eq. (8)) 


S A S D 
S A + S D 


6(1 + 6 ) 

o 


D6 


(6 + s 0 r 


1 - s. 


d 2 D 

_X_ = 
2 


( 8 ) 


Equation (8) is solved iteratively for 6. Once 6 is determined, G^ can be 

calculated from equation (1), (2), (3), and (5). Note that <5 is independent 

of delamination length. Earlier it was shown that for constant 6, Gj. and 

-4 -2 ^ —4 

Gji vary as a and varies as a . Hence, G^ varies as a and 

-2 

the corresponding applied load varies as a . The corresponding value of G^ 
-4 

also varies as a . These observations will be of special interest later when 
examining the fatigue data. 

The last parameter to be examined is the ratio of flexural-to-axial 
stiffness (i.e., D/S^) of the delaminated region. The buckling load for the 
region is linearly related to D. Prior to buckling, the load in the delami- 
nated region is linearly related to Sp. Hence, for thick specimens with a 
thin delaminated region the applied load which causes buckling depends on the 
ratio D/Sp. Delaminated regions are less prone to buckle if they possess low 
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axial stiffness and high flexural stiffness. But for homogeneous materials or 
unidirectional orthotropic laminates, this ratio is simply 



Hence, the applied load which causes buckling is independent of the material 
properties of the buckled region. However, for multi-directional laminates 
the ratio D/S d depends on both the lamina properties and the stacking 
sequence. For example, the value of S D for a [O 3 ] laminates is approximately 
1.4 times that for a [0/90/0] laminate, but the value of D is essentially 
the same for both laminates. Consequently, for a thick laminate with a 

delamination, buckling occurs at a lower applied load if the delaminated region 
is [ 0 3 ] rather than [0/90/0]. 

Comparison of Analysis and Fatigue Data 

The roles of Gj and in delamination growth were investigated by 

comparing calculated values of Gj and G n with measured growth rates. 

Fatigue data for six specimens from reference 3 were used for comparison with 

analytical results. Three of the specimens were tested at (P ) = 33 KN 

^ max 

and another three at (P^ = 30 KN. The results are presented in figure 15. 
da , 

dN decreased rapidly with delamination growth. Both curves are approximately 
linear with a slope of - 4 , hence 

da -4 

dN La (10) 

Figure 16 shows — vs. the maximum values of G ] . and G I]; during 
fatigue cycling at a maximum compressive load of 33 KN. Note that the growth 
rate is largest when Gj is relatively large. But slow growth is observed 
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even when Gj is very small. The mode II strain energy release rate 
Gji changes little after initial delamination growth. Since G 
remains large, delamination growth likely is driven by G . 

Two delamination growth criteria were examined to determine whether they 

could predict the observed growth rates. The first growth criterion examined is 
given by 


da 

dN 



(ID 


where Z and n are constants. Eqn. (11) is evaluated at the point in the 
load cycle when G^ is maximum. The load range and delamination lengths are 
such that for almost the entire test, the maximum G during each load cycle 

A **■ 

is Gj, which was obtained by solving equation (7). Recognizing that G 

A ^ 

decreases as a , as shown earlier, equations (10) and (11) can be solved for 

n; the result is n = 1. This is in strong contrast to published values of 

n 2 15-20 for double cantilever beam fatigue tests (ref. 3). Apparently, 

equation (11) is not a valid growth criterion for the specimens considered. 

Next a growth criterion was considered which includes both G and G 

1 II* 

If we assume there is no synergistic interaction of Gj and G (i.e., the 
effects are separable) , then 

dN = f l ( V + f 2 (G H J (12) 


where and are functions of Gj. and G respectively. 

From the double cantilever beam data just discussed, we know that f 
is extremely sensitive to Gj (i.e., ^-G* 5 ). Since Gj decreases 
rapidly with increasing "a", f^ must also decrease extremely fast as "a" 
increases. In fact, f^ would not contribute noticeably to da/dN after 
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the initial growth. Hence delamination growth appears to be driven by G n 
alone. Accordingly, it was assumed that the growth criterion should be 
evaluated when is maximum, i.e., at peak load. However, earlier it 

was shown that for long delaminations the crack tip closes and produces com- 
pressive a y stresses when the cyclic load is maximum. The compressive 
stress probably reduces the effect of G^ on delamination growth, but it 
was not clear how to account for this stress. Two approaches were tried: 

(1) Ignore the compressive normal stress and set f = 0 when the tip 
closes or (2) let take on negative values after the crack tip closes. 

If we choose to set ^ = 0 when the crack tip closes, then 

dN = W (13) 

for virtually the entire test. Figure 14 showed that first increased 

then decreased slightly as the delamination extends. In the experiments two 
load levels were used: (P^ = 33 KN and 30 KN. Figure 14 shows that for 
2a > 25 mm, the minimum value of G^ for (P T ) max = 33 KN is greater than 
the maximum value of G n for (P^ = 30 KN. Hence, equation (13) would 
predict that for 2a > 25 mm, the minimum da/dN for the higher load should 
exceed the maximum da/dN at the lower load. Figure 16 shows this is not the 
case. Hence, equation (13) is not valid. 

If we select a function that becomes negative when the crack tip 

closes, then we (analytically) allow compressive normal stresses at the crack 
tip to retard delamination growth due to G n . Since the compressive crack tip 
stresses increase as the delamination grows, such a function ^ would predict 
a decrease in growth rate with increased delamination length. Although this 
prediction agrees with the data trend in figure 15, more tests are needed to 
verif y or disprove this interpretation. 
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Despite the complexity of the growth behavior, two trends were clear: 

(1) high growth rates were observed only when Gj- was large, and (2) slow 
growth was observed even when G x was negligibly small; apparently, Gj 
alone can drive delamination growth. 

CONCLUSIONS 

Analysis and experiments were used to study instability-related delamina- 
tion growth in a fatigue specimen with a through-width delamination. To per- 
form the analysis an approximate superposition analysis was developed. The 
analysis expresses Gj and G n in closed form, which can be used easily to 
determine the effects of various parameters. The analysis agreed very well 
with more rigorous solutions. 

The response of the delaminated laminate to applied loads was found to be 
very complex. Key observations are listed below. 

(1) Gjj is generally much larger than Gj. 

(2) Gj and G n usually reach their peak magnitudes at different points 

in a fatigue cycle. G n always reaches its peak value at maximum 
load. 

(3) High delamination growth rates were accompanied by large values of G j. 

(4) Slow growth rates were observed even when Gj was negligibly small. 

This growth apparently was due to a large value of G . 
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APPENDIX 


A strength of materials analysis is described herein for the configuration 
in figure 1. 

The configuration is subdivided into four regions, as shown in figure 2. 
Because of symmetry, only half of the laminate is modeled. The laminate is of 
width b. The following assumptions are made: 

1. Regions B and C are perfectly bonded. Regions A and D are unbonded. 

2. Regions A, B, and C have constant axial strain. Hence, the force- 
displacement relations are those for a simple rod subjected to 
axial load. 

3. Region D has zero slope at both ends. 

4. Region D has an initial sinusoidal imperfection of peak magnitude 

The initial shape is given by 


v(x) I initial 



COS 7TX 

a 




(Al) 


where v(x) - the distortion in the y direction. 

To describe the nonlinear behavior of region D, equations (A2) and (A3) for 
post-buckling of a column were used. 

t> _ it 2 D 6 

D 26+6 (ref. (A2) 

a o 

a " a " 16a ^ + 26<S o) + ~sf (ref * 5 > (A3) 


where 6, a, a, and P D are peak lateral deflection, axial length before and 
after deformation, and load, respectively. Equations (A2) and (A3) were 
derived using strength of materials analysis of a column. 
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To combine regions A, B, C , and D, equilibrium and compatibility 

conditions must be considered. The equilibrium condition for the av-fal force 
is 


P. + 


P D' 


P*> + P, 


(A4) 


Compatibility requires the shortening of regions A and D to be identical. 
Hence , 



(A5) 


Equations (A2) to (A5) can be combined to obtain the governing equation for 
the laminate. 


7T 2 S, 


P T = 


U 6 


< 6 2 + 


266 ) 
o 


d < s a +s d> 

S A S D 


6 + 6 


(A6) 


For a specified load P T , equation (A6) is solved using a Newton-Raphson 
technique to obtain 5. can then be calculated using equation (A2) , From 

static equilibrium, the moment acting on the delaminated region at the crack 
tip is 



6 o> 



(A7) 


The force P^, is found from rule of mixtures as 

p = Sp P 
0 S A + S D T 


(A8) 
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Fig. 2- Nonlinear configuration (a) transformed into linear configuration (d) 
with two nonlinearly related loads, (P c - Pq) and M. 





Fig. 3— Typical fi 
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Fig. 4— Variation of axial strain through thickness of buckled region 
(2a = 25 mm), nonlinear analysis. 
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Fig. 5— Evaluation of approximate superposition analysis ((Pq - P^) and M 
from nonlinear finite element analysis). 
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Fig. 6- Evaluation of approximate superposition analysis ((P c - P D ) and M 
from strength of materials analysis). 
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Fig. 7- Compaction of results from geometrically nonlinear analysis. (Curve 
is visual fit of data.) 
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Fig. 9- Effect of (P c - P D ) and M on Gjj (2a = 25 mm). 
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Fig. 13- Effect of delamin 
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Fig. 16- Comparison of Gj, Gjj, and delamination growth rate for maximum 
fatigue load of 33 kN. 
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